Implementation of multipartite unitary operations with limited resources 
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A general method for implementing weakly entangling multipartite unitary operations using a 
small amount of entanglement and classical communication is presented. For the simple Hamilto- 
nian ® CTz this method requires less entanglement than previously known methods. In addition, 
compression of multiple operations is applied to reduce the average communication required. 
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Much of quantum information processing relies on the 
ability to perform operations between subsystems. In the 
case where these subsystems do not directly interact, it is 
necessary to achieve the operation by indirect means. For 
example, one could use a third subsystem which interacts 
with both, or use entanglement. In the case where the 
operation is only weakly entangling, it would be useful 
to be able to implement it using only a small amount of 
entanglement. 

A SWAP operation may be applied using two tele- 
portations P], requiring two ebits of entanglement and 
two bits of communication in each direction. Similarly, 
a CNOT may be applied using one ebit and one bit of 
communication in each direction 0. For more general 
two-qubit unitaries, weakly entangling unitaries may be 
implemented using a small amount of entanglement 
Ref. 3j gives a scheme to implement operations of the 
form e~*"'^="®'^=" using average entanglement of 5.9793q!. 
This may also be applied to more general two-qubit uni- 
taries. Alternative schemes for operations of this form, 
but allowing a significant probability of failure are pre- 
sented in Refs. 0^0- A scheme for Hamiltonians which 
are a tensor product of self-inverse operators was given 
in Ref. 0. 

Here we present a general method of implementing 
weakly entangling multipartite unitary operations using 
a small amount of entanglement or communication. The 
situation considered is that there are N parties, each with 
a number of subsystems, and operations performed by 
the parties are assumed to have negligible cost. The uni- 
tary we wish to implement acts upon one subsystem in 
possession of each of the parties, and we quantify the 
entanglement and classical communication between the 
parties. 

This paper proceeds as follows. The general problem 
of implementing operations using entanglement and com- 
munication is presented in Sec. ^ Then an efficient 
method for implementing evolution under multipartite 
Hamiltonians of the form cr®^ is presented in Sec. lIIII A 
method of using compression over multiple operations to 
reduce the communication required is given in Sec. IIVI 
A method of generalising these results to tensor product 
Hamiltonians is given in Sec. and further generalisa- 
tions are given in Sec. IVII Conclusions are given in Sec. 

IviTl 



The most general method of implementing an oper- 
ation using entanglement and classical communication 
may be described as follows. The N parties share an en- 
tangled resource state, and perform many rounds of the 
following process. Each party performs a local unitary 
followed by a partial measurement. The parties then 
communicate classical information about the measure- 
ment results, and repeat the process with local unitaries 
and measurements possibly depending on the transmit- 
ted information. 

In general we consider many implementations of U, 
and regard the scheme as successful if the minimum fi- 
delity of the output state approaches one as the number 
of copies of U approaches infinity. We may define a suffi- 
cient rate vector i? as a vector of measures of multipartite 
entanglement and classical communication such that, for 
all e > there exists an n such that n implementations 
of U may be performed using average entanglement and 
classical communication no larger than the components 
of R, and minimum fidelity at least 1 — e. 

The scheme of Cirac, Diir, Kraus and Lewenstein 
(CDKL) g is of this form. Joint measurements between 
an entangled resource state and the target system are 
made. The correct measurement result means that the 
operation is correctly performed. If failure occurs, an- 
other resource state is used, and measurement success 
corrects for the previous failure and correctly implements 
the unitary. For repeated failure this process is repeated 
until the operation which is required may be performed 
locally. 

In the following we improve on the CDKL scheme in 
three ways: 

1. We reduce the entanglement consumption. 

2. We reduce the communication in most directions to 
be the same as the entanglement consumption. 

3. We apply the scheme to more general Hamiltonians. 



III. ENTANGLEMENT CONSUMPTION 

First we give an improved scheme which reduces the 
entanglement consumed. We consider the Hamiltonian 
af^ for N parties, so the unitary we wish to apply is 
U{a) — exp{iacrf'^). This is the multipartite generali- 
sation of the Hamiltonian considered by CDKL. (We use 
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here rather than <Jx, but these are equivalent under 
local unitaries.) 

The N parties use the resource state 



|V(/3)) =cos(/3)|0) 



■isin(/3)|l) 



(1) 



then apply a four step process. In the following we use 
the "stator" formalism, introduced in Ref. 0- A stater 
is a hybrid state-operator object, and acts upon a state 
by applying its operator component to the state, and ap- 
pending its state component. That is, the stator \4>) ®U 
acts on as 



(2) 



Linearity is used to obtain the action of stators which are 
sums of terms of the form 10) ® J7. The four steps are as 
follows: 

Step 1. Each party applies a controUed-Z operation on 
the target system. This yields the stator 



cos(/3)|0)^'" ® 1^'" +isin(/3)|l) 



(3) 



Step 2. Parties 1 to — 1 apply Hadamard operations 
followed by computational basis measurements on their 
component of the resource state. The resulting stator is 



cos(/3)|0) «) 1^''' ibisin 



1) 



(4) 



where the sign is -I- if the number of measurement results 
equal to 1 is even, and — otherwise. 
Step 3. Party N performs a correction on their com- 
ponent of the resource state based on the measurement 
results. If the number of measurement results equal to 1 
is odd, then they perform a operation on the remain- 
ing portion of the resource state. We then obtain the 
stator with the plus in Eq. 

Step 4. Party N applies a projection measurement on 
their component of the resource state. Projection onto 
the state = cos(7)|0) -|-sin(7)|l) yields the unitary 
operation proportional to 



cos(/3) cos(7) (g) 1®^ + i sin(/3) sin(7) (g) af- 



(5) 



For failure, projection onto the state — sin(7)|0) — 
cos(7)|l) is obtained, giving an operation proportional to 



cos 



(/3) sin(7) (g 1®^ - I sin(/3) cos(7) (g af^. (6) 



The measurement in Step 4 is chosen such that, for the 
correct measurement result, the unitary U{a) is applied. 
This will be obtained provided 



tan(/3) tan(7) = tan(a). 



(7) 



For the incorrect measurement result, the unitary is of 
the form U{a'), but with a' ^ a. To obtain the correct 
result, we now wish to implement the unitary U{a — a'). 
Combined with the previous incorrect unitary J7(a'), this 



2p 

1.8- 

1.6- 

1.4- 

I 1.2- 
E 

- 1 
c 

a 

|0.8- 
0.6- 
0.4- 
0.2- 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.i 

a 



FIG. 1: The average entanglement consumed to implement 
the operation U{a) as a function of a. The case where the 
resource states have been numerically optimised is shown as 
the solid line, and the case for the scheme of Ref. Q is shown 
as the dashed line. 



will result in the desired U{a). To achieve this, one re- 
quires a new resource state \'ip{f3')) (where (3' is not in 
general equal to (i) and repeats the process. 

If the incorrect measurement result is obtained repeat- 
edly, one can obtain the correct unitary by using the 
resource state (lO)**^ -t- \l)^'^)/\/2. One can determin- 
istically perform any unitary of the form U{a) with this 
resource state. To see this, consider steps 1 to 4 with 
(3 = 1^/2 and ^ — a. Then, in the case of failure one 
simply applies the local unitary [/(7r/2). 

In the following the terminology "stage" is used to 
mean steps 1 to 4. The notation (3i is used for the pa- 
rameter for the entangled resource state used in stage I 
(i.e. after I — 1 failures). In addition, ai is used for the 
parameter for the unitary we wish to implement in stage 
I. That is, ai = a, U{a2) is the unitary we require for 
the correction after one failure, and so forth. We also use 
L for the maximum total number of stages. 

If iV = 2, a = 7r/2?i and j3i — ai, this method is 
similar to that of CDKL. At each stage the probability 
of success is equal to 1/2, and for repeated failure we need 
to implement U{7t/2), which may be performed locally. 
The resource state is different from that used by CDKL, 
but the average entanglement consumed is identical. 

To improve upon the entanglement consumption we 
adjust the entanglement of the resource states. By using 
resource states with larger entanglement, the probability 
for success is increased, and the average entanglement 
consumed is reduced. It is straightforward to numeri- 
cally optimise for the resource states that minimise the 
average entanglement consumed for given a. The result- 
ing average entanglement is plotted in Fig.^as a function 
of a. 

The entanglement consumed for the CDKL scheme is 
also shown in Fig. ^ The CDKL scheme is explicitly 
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given for values of a of the form a = 7r/2n. However, it 
is trivial to apply it to other values of a by making an 
expansion a = X]„ an7r/2", with a„ e {0, 1}. The entan- 
glement consumed may then be determined by summing 
it over that for the individual terms. This is how the 
values shown in Fig. ^ were determined. 

For all values of a (except the trivial a = 0) the entan- 
glement consumed for the numerically optimised scheme 
is less than that for the CDKL scheme. For small a 
the entanglement consumed is approximately 5.6418a, 
as compared to 5.9793a for the CDKL scheme. 

To show that the entanglement consumed is approxi- 
mately 5.6418a in the limit of small a, we may use the 
following approach. First we determine the entanglement 
required to implement C/(A) for values of A over a range 
of a factor of two. Here we take A g [7r/2^°, Now 
for values of a below 7r/2^°, we select an A in the range 
[7r/220,7r/2i9), such that a = A/2". Then, rather than 
numerically optimising for the best intermediate entan- 
gled states, we select the intermediate entangled states 
with f3i = a2'~^ for 1 < / < n, so (3i = ai. If a suc- 
cessful measurement has not been obtained before stage 
n + 1, then the correction required is U{A); this may be 
achieved using the numerically optimised scheme. 

The average entanglement consumed is 

n 

£(a) =£(A)2-"+^2i-'ii;(aO, (8) 
(=1 

where we use £{a) for the entanglement consumed to 
implement U{a), and E{a) for the entanglement of the 
resource state \'tp{a)). We have E{a) = /i[sin^(a)] where 
h is the binary entropy function h(p) = —p\og2{p) — (1 — 

P) log2(l -P)- 

This entanglement measure is the entropy of the re- 
duced density operator for one subsystem. This is a con- 
sistent entanglement measure for Schmidt decomposable 
multipartite pure states. It is possible to apply similar 
methods to the bipartite case to perform entanglement 
concentration and dilution between these states and the 
standard states (|0)®^ + |1)®^)/V2 0. 

The first term in Eq. (jS]) is the probability for n fail- 
ures (2~") multiplied by the entanglement required to 
implement U{A). The sum is the entanglement of the 
state with /? = aj for steps 1 to n multiplied by the 
probability. Eq. ((SJ gives the upper bound 

£(a) ^ g(A)+Er=,2^i^(A2-^) ^ 
a ~ A 

The expression on the right-hand side (RHS) of Eq. 
(PI is independent of n, and only depends on A. This 
expression is plotted as a function of A in Fig. 13 (the 
irregularity appears to be due to finite precision) . It can 
be seen that this expression does not exceed 5.6418 for 
this range of A. Thus we find that, in the limit of small 
a, the entanglement consumed need not exceed 5.6418a. 

This proof is not fully rigorous because it is not possi- 
ble to determine the RHS of Eq. for all values of A 
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FIG. 2: Upper limit on £{a)/a [the RHS of Eq. ©] minus 
5.6418 as a function of A. 



in this range. For rigour, note that a may be expanded 
as a = J2n °'ri7r/2", and the simulation of U{a) may be 
achieved by simulation of the appropriate [/(7r/2"). In 
practice, this method is more complicated, however. 

Another issue is that the above analysis is in terms 
of the expectation value of the entanglement consumed, 
given that one is provided with entangled states on de- 
mand. Alternatively we may consider many implemen- 
tations of the unitary with a fixed set of entangled re- 
source states. This is consistent with the definition given 
in SecHH 

Let p{l) be the probability of requiring the entangled 
resource state in stage I (i.e. the probability of / — 1 fail- 
ures). Now, for all e,d > 0, there exists an M such 
that the probability of the number of entangled resource 
states required in stage I exceeding Mp{l){l + 5) is no 
more than e/L. Therefore, let us consider M implemen- 
tations of U{a) with [Mp{l){l + S)\ copies of resource 
state \tp{/3i)) at each stage. Then the total entanglement 
does not exceed M£{a){l + S). 

In the case that the number of resource states exceeds 
that provided in one of the stages, let the output state be 
Pfaii- Otherwise the correct output state |0) is obtained. 
After discarding the measurement results, the average 
output state is of the form 

P = PfailPfail + (1 - P{m\)\4'){(t>\- (10) 

Because the total probability of requiring more than the 
number of resource states provided in stage I is no more 
than e/L, the total probability of failure does not exceed 
e. Therefore the fidelity must be at least 1 — e. Thus 
we find that, given sufficiently large M, the average en- 
tanglement consumed is arbitrarily close to £{a) with 
fidelity arbitrarily close to 1. 
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IV. CLASSICAL COMMUNICATION 

The next issue to consider is that of the classical com- 
munication required. At stage I the set of entangled re- 
source states is of the form 

I*,) ^ [cos(A)|0)®^ + zsin(A)|l)n^*''', (11) 

where Mi is the number of entangled states used in stage 
I, and is equal to [Mp{l){l + S)\. We may alternatively 
express this state as 

l*0=E/^'|i)' (12) 

i 

where i = (zi, . . . , jm,), A^i = Om sin'"* (A) cos^"'™ (/?;) 
and |i) = \ii) (g) . . . (Ki lu/;)- 

Now we make an approximation to the state j'l';) by 
retaining only typical sequences of the i. Denoting the 
set of typical ihy Si, the approximate state is 

|^z)(xE/.i|i)«^. (13) 
ies, 

It is a standard result for typical sequences that, for all 
e,S > 0, there exists an M; such that the fidelity is at 
least 1 — e/L with the number of elements in Si no more 
than 2^^'^(A)(i+*). 

Because fidelity does not decrease under completely- 
positive trace-preserving (CPTP) maps, the fidelity of 
the final state with the output for the exact state must 
be at least 1 — e/L. In the case where we use approximate 
resource states with typical sequences at each stage, the 
final output must have fidelity at least 1 — e with the state 
obtained with the exact resource states, p. The fidelity 
with the desired state, \(j>), is then at least 1 — 2e. 

Now, to reduce the average communication required, 
we replace Step 2 with a joint measurement in the Fourier 
transform basis on the states |i) for typical sequences 
by parties 1 to iV — 1. The communication required for 
the measurement result does not exceed MiE{/3i){l + S) 
bits. Using the expression for Mi, this does not exceed 
Mp{l)E{pi){l + S)^. 

There are a number of different ways this communi- 
cation may be performed. The communication may be 
performed from party 1 to 2, then the sum (modulo H^;!!) 
communicated to party 3, and so forth up to the com- 
munication to party N. Alternatively the communication 
may be performed directly from each party (1 to A'^ — 1) 
to party N . Many other combinations of communication 
are possible, provided it is possible for party N to deter- 
mine the sum (modulo \\Si\\) of the measurement results. 

After this communication is performed, we simply per- 
form Step 4 as above. Now it is necessary to communicate 
the positions of the measurement successes from party N 
to the other parties. It does not appear to be possible 
to reduce this communication via the typical sequences 
approach as it does for the other communication. 

The average total communication from each party 1 to 
A^-l does not exceed J2i P{l)E{fii){l+5f = E{a){l+5f. 
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FIG. 3: Results for schemes optimised for minimum classical 
communication from the last party. This communication is 
shown as the dashed line, the corresponding entanglement 
required is shown as the dotted line. The entanglement for 
the entanglement- optimised scheme is shown as the solid line 
for comparison. 



Thus, in the limit of small 5, this communication is the 
same as the entanglement consumed. The communica- 
tion from party N will, in general, be much larger, and 
does not scale down with a. For example, in Fig. ^ the 
communication required is over 9 bits for the smallest 
(nonzero) value of a plotted. 

It is also possible to select the entangled resources 
to minimise the average communication required from 
the last party. In this case numerical results indicate 
that the optimal scheme uses an entangled resource with 
tan(/3) = tan^(Q;). For small a, the probability of failure 
is small. In the case of failure we simply implement the 
scheme with 1 ebit of entanglement and 1 bit of commu- 
nication in each direction. The corresponding classical 
communication is shown in Fig. |51 

The classical communication is still larger than the en- 
tanglement for the scheme optimised for minimum entan- 
glement consumption. However, it does go to zero in the 
limit of small a. Unfortunately the ratio of the classi- 
cal communication to a is not bounded for small a (see 
Fig. The ratio increases approximately as |loga|. 
This means that a scheme based on implementing a se- 
ries of operations of the form U{a) will still require a 
large amount of communication from the last party. 



V. GENERAL TENSOR PRODUCT 
HAMILTONIANS 

In the previous sections we have shown how to imple- 
ment evolution under a Hamiltonian of the form crf^. 
Now we show how to implement a general tensor prod- 
uct Hamiltonian of the form H = Hi ^ . . . ^ //jv- To do 
this, we apply a method similar to that in Ref. fSj. First 
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FIG. 4: The ratio of the communication from the final party 
to the interaction strength a. 



note that it is possible to diagonalise the Hj via local 
unitaries. Taking A = the diagonalised form of H 



Ha 



|diag(aij, a2j,a3,^ 



(14) 

where aij G [—1, 1], and dj is the dimension of the sub- 
system which Hj acts upon. In the following we consider 
simulation of this diagonalised form. 

This Hamiltonian may be simulated by the Hamilto- 
nian af^^ in a similar way as in Ref. 0. The chain of 
simulations used is 



H'i 



A2B2 



Hdiag/ A, 



(15) 



be achieved by simply restricting to the appropriate sub- 
space. 

Ref. 1^ gave a similar method for the case of bipartite 
tensor product Hamiltonians (though the multipartite 
case was mentioned briefly in the discussion). In the bi- 
partite case of Ref. the simulation is reversible. This 
is because the diagonal Hamiltonian may be expanded as 
a sum of local Hamiltonians and a Hamiltonian for which 
the maximum and minimum eigenvalues have the same 
magnitude. The local Hamiltonians may be ignored, be- 
cause they can be implemented locally. One can then use 
i/diag/A to simulate trf ^ by restricting to the subspace 
for the maximum and minimum eigenvalues. In the mul- 
tipartite case we can not use a similar method, because 
the additional terms will be multipartite Hamiltonians 
on iV — 1 subsystems, rather than local Hamiltonians. 

The case of self-inverse Hamiltonian evolution, as in 
Ref. 0, is particularly simple. All the aij are ±1, so it 
is possible to achieve the entire simulation chain without 
intermediate unitaries simply by restricting to the ap- 
propriate subspaces. Thus self-inverse Hamiltonian evo- 
lution is equivalent to evolution under crf^. 

In order to simulate evolution under H for a time in- 
terval t, we divide the time interval into m intervals of 
length t/m. For m sufficiently large, the chain of simu- 
lations above will be accurate. The expectation value 
for the consumed entanglement is then no more than 
5.6418i||77|| . In fact, the expectation value is approxi- 
mately equal to this, because each simulation of U{a) is 
for a <C 1. As above, we may also obtain an asymptotic 
average entanglement consumed equal to this value in the 
limit of a large number of implementations of H. Also 
the asymptotic average communication from each of the 
parties 1 ... iV — 1 may be made equal to this value. 



where H'^,^, = diag(ai j , -oi j , . . . , ad^j,-ad^j). We 
use the Aj subscripts to indicate the subsystems upon 
which we wish to implement H , and Bj to indicate ancilla 
subsystems. 

In order to perform this chain of simulations, we need 
to, in general, perform the simulation 



^®(JV-i+i) 



' ^A,^B„ 



(16) 



Note that it is necessary to consider a large number 
of implementations of the Hamiltonian for small time in- 
tervals, because in addition to the division of t into m 
intervals in order to make the simulations in Eq. I|15() 
accurate, the simulations in Eq. (|15|l require further sub- 
division of the time. However, for a given fidelity 1 — e 
required, the number of subintervals is fixed, so there is 
no problem in taking the appropriate limits to obtain the 
average entanglement and communication. 



To do this, we first append the dj dimensional an- 
cilla Bj so the j'th term in the tensor product is 
diag(l, —1, 1, —1, . . . , 1, —1). Now we take pi = {aij + 
l)/2 G [0, 1], and define the local unitaries Ui which ex- 
change the {21 — l)-th and (2^)-th basis vectors of AjBj. 

To simulate {g)^^i H'^^b^ ® cr?^^""^'^ for small time 5t, 
we simply apply Ui at time pi5t and again at time 8t. 
In this way, we may apply each successive simulation in 
the chain H15|) . The final simulation in the chain may 



We can also use the above methods to simulate more 
general Hamiltonians of the form H = X^/c-^fe' where 
each Hk is a tensor product Hamiltonian. We simply 
implement (H^e"'*^''/'")'". For large to, this simulation 
may be made arbitrarily accurate. The average entangle- 
ment consumed is approximately 5.6418t^j. ||-fffe||. For 
completely general Hamiltonians a large number of terms 
in the sum will be required, so this method will not be 
efficient. However, it should be useful for cases with only 
a moderate number of terms required. 
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VI. FURTHER GENERALISATION 

We may generalise the above methods to general mii- 
taries of the form 

[/ = ^AfcF«®yf (17) 
k 

where the V^''^ are local unitaries, and is small except 

for the k such that V^f'''' = 1. To implement this op- 
eration, we use an entangled resource state of the form 

I'^u) ^^IJ-k\k)B^(g> ■.■(E>\k)B,^. (18) 

k 

We use Bj for the subsystems for the resource state, and 
Aj for the subsystems for the system state (which we 
wish to apply U to). 

Each party j performs the controlled operation 
J2k m'^V^^\ where V^^^ acts upon Aj. As a result, 
we have the stator 

^ fik\k)B^ ® vl:'^ ® . . . ® |fc)B„ ® vl""^ (19) 

k 

acting upon the system. 

Each party performs local measurements, and commu- 
nicates the results to party N, where the phase correction 
is made. We then perform a projective measurement on 
the resource state in subsystem Bn^ with success corre- 
sponding to projection onto the state i^klk). Choos- 
ing /ifcz/^ oc Afe, U is implemented in the case of success. 

In the case of failure a different unitary operation is 
performed. This case may be corrected for by imple- 
menting the correction in the same way, or alternatively 
by teleporting the portions of the state in subsystems Ai 
to Apf-i to party N, performing the appropriate unitary, 
then teleporting these portions of the state back to their 
respective subsystems. 

For example, we may use this approach with two-qubit 
unitaries. Two-qubit unitaries may be simplified to the 
form Ulllll 

3 

C/ = ^AfcCrfc(^crfc, (20) 
fc=0 

where we take the convention that cto = 1- In the case 
of the incorrect measurement result, the error operation 
is again of this form, and we wish to implement an op- 
eration of this form for the correction. Repeating the 
process, we may implement the unitary in a more di- 
rect way than by implementing a sequence of unitaries of 



the form U{a). However, calculations indicate that this 
method consumes more entanglement. 

One use of this approach would be to implement the 
operation with small classical communication in all di- 
rections. As discussed at the end of Sec. lIVI it is possible 
to optimise for minimum communication from the last 
party, however the ratio of the classical communication 
to a is not bounded. This means that, if we use this 
method for more general tensor product Hamiltonians, 
the average communication from the last party is still go- 
ing to be large. On the other hand, using the approach 
given in this section, we can ensure that the average com- 
munication required approaches zero in the limit as U 
approaches the identity. By selecting fXk oc oc \/Afc, 
we find that the probability of failure approaches zero as 
U approaches the identity. 



VII. CONCLUSION 

We have presented a scheme for implementing evolu- 
tion under general multipartite tensor product Hamilto- 
nians. This scheme requires multipartite entanglement 
of approximately 5.6418i||i/|| . In contrast, the scheme 
of CDKL Q, which applies to the O ctx interaction, 
requires entanglement of 5.9793i||-ff |j . This scheme may 
be applied both to the case of a single implementation, 
in which case this is the expectation value of the entan- 
glement consumption, or to multiple implementations. 

In the case of multiple implementations, the entangle- 
ment consumption is the average over the implementa- 
tions, given that the implementations are achieved with 
arbitrarily high fidelity. In addition, for multiple imple- 
mentations the average communication in some of the di- 
rections may be reduced to be equal to the entanglement 
consumption. However, a large amount of communica- 
tion from one of the subsystems is required. 

This scheme may also be applied to more general 
Hamiltonians. A general Hamiltonian may be expressed 
as a sum of tensor product Hamiltonians, and therefore 
may be implemented using the above approach together 
with a Trotter expansion. The drawback to this approach 
is that it will be inefRcient if a large number of terms is 
required. It is also possible to apply this approach to 
CPTP maps which are close to the identity. By adding 
an ancilla, the map may be replaced with a unitary which 
is close to the identity, which may be implemented in the 
above way. Again this method will be inefficient if a large 
number of terms are required. 



[1] C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. (1993). 
Peres, and W. Wootters, Phys. Rev. Lett. 70, 1895 



7 



[2] J. Eisert, K. Jacobs, P. Papadopoulos, and M. B. Plenio, 

Phys. Rev. A 62, 052317 (2000). 
[3] J. I. Cirac, W. Diir, B. Kraus, and M. Lewenstein, Phys. 

Rev. Lett. 86, 544 (2001). 
[4] B. Groisman and B. Reznik, Phys. Rev. A 71, 032322 

(2005). 

[5] L. Chen and Y.-X. Chen, Phys. Rev. A 71, 054302 
(2005). 

[6] H.-S. Zeng, Y.-G. Shan, J.-J. Nie, and L.-M. Kuang, 

quant-ph/0508054 (2005). 
[7] B. Reznik, Y. Aharonov, and B. Groisman, Phys. Rev. 



A 65, 032312 (2002). 
[8] C. H. Bennett, S. Popescu, D. Rohrlich, J. A. Smolin, 

and A. V. Thapliyal, Phys. Rev. A 63, 012307 (2001). 
[9] A. M. Childs, D. W. Leung, and G. Vidal, IEEE Trans. 

Inf. Theory 50, 1189 (2004). 
[10] N. Khaneja, R. Brockett, and S. J. Glaser, Phys. Rev. A 

63, 032308 (2001). 
[11] M. A. Nielsen, C. M. Dawson, J. L. Dodd, A. Gilchrist, D. 

Mortimer, T. J. Osborne, M. J. Bremner, A. W. Harrow, 

and Andrew Hines, Phys. Rev. A 67, 052301 (2003). 



